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^ Abstract 

D We review basic theoretical concepts to describe nuclear collisions at ul- 

trarelativistic energies. We discuss relativistic strings, Gribov-Regge theory 
\^ (GRT) of hadronic interactions and string models based on GRT, and finally 

generalizations to nucleus-nucleus scattering. 
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^ 1 Introduction 

The early universe was probably a hot and dense "fireball" of quarks and gluons, 
before, due to expansion and cooling, hadrons emerged. Presently, there are con- 
siderable efforts to create such a "quark gluon plasma" (QGP) in nucleus-nucleus 
scattering at ultrarelativistic energies 1 GeV per nucleon) [1]. 

In ultrarelativistiv collisions there is certainly enough energy available to produce 
high enough energy densities for a QGP, provided all the energy is used to heat up the 
system. This is, however, not the case. We know that the nuclei are to some extent 
transparent, they go through each other by keeping a large fraction of their original 
momentum. But still, the nucleons do loose energy, which shows up as baryon-poor 
matter in the central region. To be more precise: the system has roughly the form of 
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a cylinder which expands essentially longitudinally, the forward and backward front 
moving almost with the velocity of light. The forward and backward region of the 
cylinder are baryon-rich, the central region is baryon-poor. On the other hand, the 
energy density is largest in the central region, and the big question is whether the 
energy density in this region is large enough to form a plasma. 

Theoretically, such questions related to the space-time structure can be invest- 
igated by using string models [2]-[6]. These models provide a full description of 
ultrarelativistic nucleus-nucleus collisions, starting really with two incident nuclei 
and not at some vaguely known intermediate stage as many hydrodynamical mod- 
els. Basic features of high energy hadronic collisions as tranparency, longitudinal 
structure . . . can be easily understood in terms of string models, even more, these 
are basic properties of string models and not just a specific choice of parameters. 
String models are also very useful to study signals of the QGP, since practically all 
signals are strongly affected by the space-time evolution of the system. 

There is quite a number of microscopic dynamical models, a fact which is at first 
sight somewhat disturbing. Looking closer, one finds, however, similarities. There 
is a whole class of models (VENUS [2], the dual parton model (DPM) [3, 4], the 
quark gluon string model (QGSM) [5, 6]) with all the models being strictly based 
on Gribov-Regge theory and Veneziano's cylinder hypothesis (these models are re- 
ferred to as Gribov-Regge models, GRM's). This general framework is referred to 
as "the string model" in this article. The different models of this class (VENUS, 
DPM, QGSM) differ in details concerning the precise formulation of string form- 
ation and decay. A nice feature of the Gribov-Regge approach is the consistent 
formulation within relativistic quantum theory. The classical string is only used as 
a phenomenological parametrization of particle production, this is not a classical 
model. 

There is no real alternative to the string model discussed in this article. There 
is a "classical string model" [7], which seems to be quite different but nevertheless 
successful. However, also here the hadronic interaction amounts to string formation 
and decay. Lacking theoretical guidance one simly assumes something about string 
formation, which is actually very similar to the "colour exchange mechanism" of the 
GRM's. Another successful approach is the RQMD model [8], which uses classical 
trajectories and measured hadron-hadron cross sections. But also here, for energetic 
hadron-hadron interactions, a string approach is applied. Eurthermore, at high 
energies, the Gribov cross sections are used to introduce multiple scattering. So all 
successful models use at least elements of GRT, even when they are not formulated 
within this framework. 

In this article we treat only the class of models based on Gribov-Regge theory 
(VENUS, DPM, QGSM), the general framework referred to as "the string model". 
We discuss the theoretical concepts as relativistic strings and Gribov-Regge theory, 
and demonstrate how both are linked to provide "the string model of hadronic 
interactions". We also discuss the generalization to nucleus-nucleus scattering. Many 
details and references missing in this article can be found in a recent review article 
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[2]. 

2 String Dynamics 

The dynamics and fragmentation of relativistic strings are crucial ingredients of the 
model for hadronic interactions to be introduced later. We shall discuss strings in a 
general fashion in order to demonstrate that dynamics and fragmentation are almost 
fixed from symmetry requirements. This is the reason that string fragmentation 
models are much less arbitrary than one might think. 

We review classical string theory [9, 10, 11]. We discuss how to obtain a string 
action, we derive equations of motion for the space-time evolution of strings, and we 
discuss the general solution as well as conservation laws. We then treat the simplest 
possible string solution, the so called "yo-yo" string. 



2.1 A Gauge Invariant String Action 

We discuss in this subsection how to obtain a string action from invariance require- 
ments. 

A classical string is a two-dimensional surface in the four-dimensional Minkowski 
space, 

X = a;(r, cr) , (1) 

with a spacelike parameter a and a timelike one r. Of course this is only one of 
infinitely many parametrizations of this surface. A transformation 

T 

cr 

from one parameter space to another is called a gauge transformation^ and the 
group of such transformations is called a gauge group. One assumes that the string 
action should not depend on the parametrization, so gauge invariance is a necessary 
requirement. Further restrictions should be locality and covariance. Concerning 
the question of gauge invariance, it is useful to relate a metric g to a certain string 
parametrization via 

gap = daX^dpXf, , (3) 

where a and f3 assume the values 1 and 2, and where we used di = and 82 = -§^- 
By taking "dot" and "prime" as abbreviations for ^ and the metric can be 
written as 



ry* ^ ry* ry* ^ ry* ^ 

fAy fAy tAy tAy 



(4) 



How does the metric g transform under the gauge transformations given in eq. (2)? 
Defining the two component variable ^ via 



^1 = r , 6 = cr , 



(5) 
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and using 
we get 

gafj = dax^dpx^ (7) 

= d.x^'dc^i.djx^dpij (8) 

= dalfg.jdpij . (9) 

Since the components of the Jacobi matrix M of the gauge transformation eq. (2) 
are given as 

Ma, = dkia , (10) 

we can write eq. (9) in matrix notation as 

g = M^gM . (11) 

This leads to the identity 



det^l = ^|det^| |detM| . (12) 
On the other hand we have 

d'^i = I det M\ (P( , (13) 



which together with eq. (12) immediately suggests that a !^ integration over ^| det g\ 
is invariant under gauge transformations, 



/ = j ^\Aeig\<f( = j ^\Aeig\<f( = I . (14) 
Writing the integral / explicitly as 

1=1 \J{x'xY — x'^i^ drda (15) 



shows that / is also local and covariant. In fact, / is the simplest local, covariant, 
and gauge invariant expression, and so / is a very attractive candidate for a string 
action. Therefore we define the action of a relativistic string to be [9] 



S = I L drda , (16) 

with 



L = —K \J — det g = —k y (x'i)^ — x'^i^, (17) 

where we used | det g\ = — det g. We will see later that the proportionality constant 
n can be identified with the "string tension", the energy per unit length of the string. 
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2.2 Equations of Motion, Conservation Laws 

We rewrite the action defined in the last subsection more explicitly as 



S = r dr r da L , (18) 
Jti Jo 



wi 



ith 



L = —K \l{x'xY — x'^i^ . (19) 



We use the convention CTmin = and CTmax = tt, where tt is an arbitrary number at 
the moment. The symbols ri and T2 represent initial and final times. To obtain the 
equations of motion, we require 

5S = {) (20) 

under infinitesimal variations 8x{(t,t) of the string surface. We find the equations 
of motion 

+ (21) 
dr dxf^ da dx'^ ' 

and the boundary conditions 

dL 



dx'^ 



at (7 = 0,^. (22) 



From the invariance of the action under translations one obtains conservation 
laws for energy and momentum. Defining the energy-momentum currents as 

we may introduce the string momentum in various ways, for example as 



(string) := / da , (24) 

JCr 

where integration at constant r is implied. From 8S = 0, one concludes that the 
string momentum is conserved. 

Using the currents defined in eq. (23), we may rewrite the equations of motion 
eq. (21) as 

|p;' + |p,- = 0, (25) 

and the boundary condition eq. (22) reads 

P^^ = , at (7 = 0, ^ . (26) 



Our next aim will be to solve these equations of motion. 
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2.3 Solutions of the String Equations 

To solve the equations of motion, we choose a gauge which simphfies the equations 
of motion. The orthonormal gauge 

x'x = , + x'"^ = (27) 

does so. The currents eq. (23) are now 

Pr = KX , P„ = -Kx' , (28) 

the equations of motion eq. (25) are simply wave equations 

- < = , (29) 

and the boundary conditions eq. (26) are 

x'{t,0) = x'{t,TT) = . (30) 

To completely specify the gauge, we set 

Xo = t = T . (31) 

This implies that r and a have length dimensions. For the following we use this lab 
frame parametrization^ and consider only the space components of x. Eq. (27) now 
reads 

f'f = 0, (f)^ + {xf = 1 . (32) 
The solution of eqs. (29, 30) is 

x{t,a) = ^[y{t + a) + y{t-a)] , (33) 

where y[t) is obviously the trajectory of one endpoint, x(t,0), called the directrix. 
The directrix has to be periodic, 

2P 

y{t + 2tt) - y{t) = — , (34) 

K 

where P is the string momentum (to be shown later). Eq. (33) has the follow- 
ing meaning: each point on the string may be obtained by a simple geometrical 
construction once the directrix y[t) is known. 



2.4 The Yo-Yo String 

A simple but important example is the so-called yo-yo string, characterized by a 
one-dimensional directrix with one period consisting of two linear segments. In case 
of a one-dimensional directrix, straight lines with a tilt of 45° against vertical (in 
space-time) are mandatory, because the string end (represented by y{t)) moves with 
the velocity of light (because of eqs. (32, 30)). Erom eq. (33) it is clear that 
the corresponding string is a simple straight line [qq in fig. 2) stretched between 
directrix and antidirectrix. 
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antidir. 



E 




Figure 1: A string with its dir- 
ectrix and antidirectrix. The dir- 
ectrix segment DAE defines the 
string piece AB and the anti- 
directrix FCG the string piece 
CB. 




3 String Fragmentation 



We discuss the rules for string breaking in the framework of classical relativistic 
string theory in general and in particular for yo-yo strings. Although, in classical 
string theory, the time evolution is fixed once a string breakpoint is known, the 
determination of locations of breakpoints requires further input. For this purpose 
we employ the same symmetry arguments which led earlier to the string action. This 
procedure leads to the "area law" [12, 13, 14], i.e. the probability dP of a string to 
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Figure 3: Breaking of a "yo-yo" string: we first determine the segments of the 
(anti)directrix corresponding to the string pieces AB and CB; these segments are 
then periodically continued into the future. 

break within a given area element d^A is given as 

dP = {1 - P) a d^A , (35) 

with the "break probabihty" a as a parameter. Knowing the breakpoint, it is clear 
know how to proceed. As for the action, we assume locality. If a break occurs at 
x[t^ cr), we have to make sure that for the future as well as the past we have periodic 
(anti-) directrices, and that the directrices for future and past match properly in 
the present. The only way to ensure this is to periodically continue the directrix 
corresponding to one string piece and the antidirectrix corresponding to the other 
string piece into the future. This fully determines the time evolution of either string 
piece also for all the future, at least untill the next break. 

Let us now discuss these "cutting rules" for a yo-yo string (see fig. 3). Without 
interaction, the string stretches between directrix (t, y{t)) and antidirectrix (t, y{t)) = 
(t, a;(t, tt)). Let the point B = [t^x[t^cr)) be a breakpoint on the string at time 
t, dividing the string into two segments AB and BC with A = (t,a;(t,0)) and 
C = (t, x[t^ tt)). The directrix and antidirectrix corresponding to these segments are 
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DEF with 



and GHI with 



Using 



D = {t-a,y{t-a)) , (36) 
F ={t + a,y{t + a)) (37) 



G =(t-(^-a),y(t-(^-a))) , (38) 
/ = (t + (^_a),y(t + (^-a))) . (39) 



y{t) = y{t-n) + - , (40) 

K 



we verify easily that, after the appropriate shift, the segments DFF and GHI 
provide a full period of the unperturbed string. As discussed earlier, we obtain 
the directrices of the two segments after the break by continuation of DFF (— > 
DFFS ■ ■ ■) and of GHI (— > GHIN ■ ■ ■ )■ The corresponding antidirectrices can be 
easily constructed from the relation 

y{t) = \{y{t + ^) + y{t-^)) (4i) 

between directrix y and antidirectrix y, and we get BKM-- - and BQT---. We 
realize the identities 

\\BJ\\ = \\JK\\ , \\HJ\\ = \\JI\\ (42) 

and 

||5P|| = ||Pg|| , ||PP|| = ||PF|| , (43) 

which provide a very simple procedure for actually constructing the new directrices 
in numerical applications. 



4 Gribov-Regge Theory (GRT) 

Gribov's multiple scattering theory of ultrarelativistic hadronic interactions, referred 
to as Gribov-Regge theory (GRT), is the theoretical basis of the string model of 
hadronic/nuclear scattering to be introduced later. 

We briefly discuss the Pomeron. We then introduce an expression for the amp- 
litude of elastic hadron-hadron scattering due to multiple Pomeron exchange. We 
calculate the elastic and, via the optical theorem, the total cross section. We discuss 
the AGK cutting rules, which provide a technique to calculate the inaginary part of 
elastic amplitudes. Finally we apply the AGK rules to expand the total cross sec- 
tion as (Ttot = Z]Cm, with "topological cross sections" (7^ referring to m elementary 
inelastic processes (m cut Pomerons). 
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1 2 3 




Figure 4: 

n Pomeron exchange. 



4.1 The Pomeron 

In the next subsection we introduce a multiple scattering theory of hadron-hadron 
interactions, implying multiple exchange of elementary objects called Pomerons. 
From general considerations of the high energy limit of elastic amplitudes, one para- 
metrizes the amplitude associated with Pomeron exchange as 

~ ?^ (44) 

with s and t being the Mandelstam variables s = [pi +^2)^ and t = [pi — ps)^, with 
Pi} P2 (P3, P4) being the momenta of the incoming (outgoing) hadrons. 

The nature of the Pomeron in terms of quarks and gluons is still not known. 
Initially, Pomerons were thought to be ladder diagrams (gluon ladders), or gluon 
networks of cylindrical topology, but QCD calculations are not conclusive. We adopt 
Venezianos picture of a Pomeron being a cylinder of gluons and quark loops (being 
a generalized gluon ladder). 

4.2 The Multi— Pomeron Amplitude 

Starting point is the following expression for the elastic amplitude: 

A2^2{s,t) = J2A.,{s,t), (45) 

n 

with 

:n-l l-n ,. n 

t) = ^ / n d\ S^'Xk - ki) N^h ■■■K) D{s, kl)--- D{s, kl) , 

nl J ^ 

(46) 

representing n Pomerons exchanges (fig. 4) . The variables k and ki represent 
transverse momenta. The Pomeron Green's function is 

/ „ \ a(-fc2)_l 

D{s,e) = r^[-J . (47) 
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Introducing a rapidity gap y via y = Ins/so and defining 

A := a(0) - 1 , (48) 

we find 

D{s, k^)^ri exp(Ay) exp [ - a'(0) y k^] , (49) 

witli exp(Ay) being 1, if we liave an ideal Pomeron witli a(0) = I. As we are going 
to argue later, the data, however, suggest that a(0) is slightly larger than I. The 
simple form of the Pomeron Green's function is a pure assumption (of Regge pole 
dominance), only justified by the success of the theory, and it is, therefore, most 
important to justify this assumption within QCD. 

Eqs. (46, 49) are the basis for applications to be discussed in the following. 



4.3 Elastic Scattering and Total Cross Section 

As a first application of the Gribov-Regge theory, we consider elastic scattering. 
The amplitude for elastic scattering is given as 



A{s,t) = Y,A^{s,t) 



(50) 



n=l 



where An represents n Pomeron exchanges, and is given by eqs. (46, 49). Assuming 
factorization of the vertex function, 



7 , 



Yl N{k 



I ) 1 



8 = 1 



we get 



A(s,t) = — / c/% exp{ikb) ^{s,b) 

4:TT 



ith 



wi 



and 



7(3,6) = i{l-exp[-Ccu(3,6)]} 

674 



A^o exp(Ay 

CU(S,Oj = ; ; exp 



i?2 + a'y 



i?2 + a'y 

Eqs. (52, 53, 54) may be used to calculate cross sections as 

(jtot = 8^ImA(s,0) = J (/%2Re7 



an 



d 



dkHn \A\ 



(51) 

(52) 
(53) 
(54) 

(55) 
(56) 



There are five free parameters, A^o, A, R^, a', and C, which are fixed by com- 
paring with data. The increase of crtot('5) and crin(s) with s as well as many elastic 
scattering data can be nicely reproduced, not so for A = 0. 
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^ 

4.4 The Abramovskii— Gribov— Kancheli Cutting Rules 

Due to the optical theorem, the discontinuity (or imaginary part) of the elastic amp- 
litude is related to inelastic processes. One therefore investigates discA in order to 
study inelastic scattering. The Abramovskii-Gribov-Kancheli (AGK) Cutting Rules 
[15] provide a technique to express discA (or ImA) in terms of discG, representing 
elementary inelastic processes associated with the exchange of a single Pomeron. 

We start with the elastic amplitude for multiple Pomeron exchange, given in eq. 
(46), which we write as 

/n 

7=1 

with Gi = G{s^kf) = N[kf)D[s^kf). The Cutkoski cutting rules state for a Feyn- 
man diagram (or a sum of graphs) having the structure eq. (57): 



If 11 

- discA„ = Y. H iGa H -G} n - discG^ = Y.^r. 

* cuts left riftht * cut * m 



(58) 



cuts left right 

of cut of cut 



where m represents the number of elementary inelastic interactions ("cut Pomeron- 
s"). One finds analytic expressions for Anm [2]. 

4.5 An Expansion of the Total Cross Section 

We are now going to use the expansion of 4- discA„ in terms of the number of cut 
Pomerons to obtain an expansion of the total cross section in terms of topological 
cross sections: atot = Y,^m- Here cr^ corresponds to the cross section of m element- 
ary inelastic processes (m cut Pomerons). We have 

1 47r 
(7tot(s) = — discr(s,0) = — discA(s,0) (59) 
2is I 

oo -j^ 

= 4^ y -discA„(s,0) , (60) 
where n refiects n Pomeron exchanges. We obtain 

oo 

^tot(.s) = XI ^-nii^) 5 (61) 
m=0 

with 

8 ^ iVo exp( Ay) r _^ '"i^ 



[l-e'^Ed '"^>0' (62) 



m z L ,„ k\ 

k=0 



which may be evaluated easily numerically. For m = 0, we get 

ao(3) = 8^iVoexp(Ay)X ^ (--) ■ (63) 



-, n n\ 

n=l 



The (Tm are strictly positive! 
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Figure 5: Cylinder diagram (gluons 
and quark loops on the back sheet are 
not drawn). 



Figure 6: Cut cylinder diagram 
(gluons and quark loops on the back 
sheet are not drawn). 



5 The VENUS String Model 

The Gribov-Regge theory provides a framework to calculate amplitudes ^2^2 for 
elastic hadron-hadron scattering, and therefore elastic and total cross sections can 
be calculated. For more detailed investigations, one needs inelastic amplitudes A^^n^ 
describing particle production. These amplitudes are, however, not calculable within 
GRT. On the other hand, we know that discontinuities (or imaginary parts) of elastic 
amplitudes are related to inelastic scattering, so one may take the expansion of 
discA2^2 (or of (Jtot) from the last section as a guideline to construct a model. This 
is exactly what is done in the VENUS model: 

• In VENUS, elastic amplitudes are calculated strictly according to Gribov- 
Regge Theory (GRT), and so are the elastic and total cross sections. 

• VENUS provides a model for calculating inelastic amplitudes, guided by the 
expansion of (Jtot in terms of topological cross sections. 

The VENUS model to be introduced in the following is closely related to the 
dual parton model (DPM) [3, 4], introduced by Capella et al., and the quark gluon 
string model (QGSM) by Kaidalov et al. [5, 6]. 

In order to formulate a model for inelastic scattering, one needs to know something 
about the nature of the Pomeron. According to Veneziano [16], a Pomeron is a cyl- 
inder, see fig. 5. Correspondingly we identify the discontinuity of the Pomeron 
propagator (or more precisely — zdiscG = 2 ImG) with a squared cut cylinder, with 
a cut cylinder shown in fig. 6. The two cutting edges of the cut cylinder are 
identified with relativistic strings, as shown in fig. 7, where the thin lines represent 
(anti)quarks and the thick lines socalled remnants, each one representing an incident 
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Figure 7: Cut cylinder corresponding to two strings. 



hadron minus the corresponding (anti)quark. So, before the interaction, (anti)quark 
and remnant of the same hadron are connected, after the interaction, (anti)quark 
and remnant from different hadrons. Therefore the term "colour exchange" is used, 
because it looks like colour being exchanged between the two (anti)quarks. 

Considering two incoming hadrons /2(p2), with f[p) refering to flavour 

(momentum), and a colour exchange between the (anti)quarks i[k) and one 
obtains two strings S+ and S-. The string 

S+:={Mp^)-z{k)\j{l)} (64) 

contains at one end the projectile remnant fi[pi) — i[k)^ being the projectile /i re- 
duced by the parton (quark or antiquark) i. This string end carries the projectile 
momentum pi reduced by the parton momentum k. The other end consists of the 
parton j with momentum /. Since this string contains the projectile remnant, which 
has a large forward momentum, we refer to this string as forward string or forward 
baryonic string (it carries baryon number f). The other string 



S. ■.= {f,{p,)-j{l) t{k)\ (65) 

contains at one end the target remnant /2(p2)— j(/), which is the target /2 reduced 
by the parton j, with the target momentum p2 reduced by the parton momentum /. 
On the other end, we have the parton i with momentum k. This string is referred 
to as backward string or backward baryonic string (it carries as well baryon number 
f ). So if, for example, /i and /2 are protons and i and j are u quarks, the two strings 
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(D) 



/i 



J 



c 
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c 
c 
c 
c 
c 



Figure 8: Nondiffractive (N) and diffractive (D) type. 



are ud-u strings, with the diquark ud being in one case the forward end and in the 
other case the backward end. 

The forward string S+ consists of a backward parton j and a forward hadron / 
minus a parton i. It might be by chance that fi — i is equal to /i + z with /i having 
survived as a singlet, or, in other words, the CE involved only a colour singlet i — i 
pair and left /i surviving as a spectator. The same arguments apply of course for 
the backward string S-. We refer to such interactions as "diffractive" (D) and to 
others as "nondiffractive" (N) (see fig. 8). We say, the hadron suffers a D-type or 
N-type interaction (we use weights w and 1 — w). 

For the full I-cylinder contribution (consisting of two strings), we have four 
combinations, NN, ND, DN, and DD, with weights (I — t/;)^, (I — w)w^ w[l — w)^ 
and t/;^, which are referred to as 



These contributions are conveniently represented as quark-line diagrams as shown 
in fig. 9. 

6 Nucleus-Nucleus Interactions 

We are going to generalize the results of the preceeding sections to nucleus-nucleus 
[A-B ) scattering in a straightforward way. 

For elastic scattering, we consider diagrams of the type shown in fig. 10, where 
the dashed line with index n represents a block of n Pomeron exchanges. The same 
nucleon may (or may not) be involved in more than one interaction. We write the 



nondiffractive scattering (NN) 

diffractive projectile excitation (ND) 

diffractive target excitation (DN) 

Pomeron — Pomeron scattering (DD) 



(66) 
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(NN) 
9- 



(ND) 
<9- 



©- 
©- 
©- 

©- 



(DN) 
<9- 



fDDl 



©- 
©- 
©- 

©- 



Figure 9: Quark line diagrams representing the four contributions to single colour 
exchange between quarks. 



elastic amplitude for A-B scattering as 

^ = E E E K''nr'' , (67) 

aif3i...ai,f3i, ni...nv 

with 

/ N^J^^-'^^" n — n ^ ^ ' (68) 

tl=\ ^1^- 'J=\ 

where represents a sum over all possible "collision sequences" ai/3i . . . 

with and being indices refering to the projectile and target nucleon involved 
in the interaction. Ordering in the sequence is irrelevant, and the double indices 
have to be pairwise different: aif^i ^ a^fSj. Eq. (68) corresponds to a sequence of u 
N-N collisions, each of them representing a "block" of several pomeron exchanges 
(see fig. f 0). Now we expand the absorptive part of i A'^^^^^'^""^" as a sum of terms 
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Figure 10: Elastic nucleus-nucleus scattering: nucleons of the projectile nucleus 
interacts with nucleons from the target nucleus via a "block" of several Pomeron 
exchanges (dashed line) for each N-N interaction. 



corresponding to out of the Pomerons being cut, 



1 



■ ni...ni^ / J n\...nv ,m\...mv 



(see fig. f f), with 



\0.\fi\...O.vfiv 



n 



/ y I ni...nv 

r T J 



lelhh 

I n 



— - n 2 ImD 



n i-^T'^^ 



n, 



n 2 ImD -2 Re ]\-D 



1 1 
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(69) 



(70) 

/in, /el, and /ib are sets of indices, corresponding to inelastic, elastic, and inter-block 
cuts: /in is the set of indices [i with > 0, and /gi, hh are partitions of /j, \ /in, 
with Ii, = {I ... I/}, which means /gi U /ib = Ii, \ /in- The sum Zl/^j/;^ is meant to sum 
over all partitions. Using the optical theorem, one obtains 
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Figure 11: A contribution to the absorptive part of the elastic amplitude for nucleus- 
nucleus scattering, related to cutting out of the "block" of exchanged Pomerons 
in the nucleon-nucleon interaction. 



with 
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(72) 
(73) 



a = l 



m :=^{l-exp[-2Cc.(6)]} 



c 

1 \2Cuj{h) 



C ml 

Performing the sum in eq. (72), one obtains 



exp - 2Cuj{b)\ . 



. AB 
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(74) 

(75) 
(76) 

(77) 



(78) 
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This formula, derived in the GRT framework without approximations (apart from 
assuming a simple form of the nucleus-nucleon vertex) is the starting point of many 
applications. Eq. (78) is often referred to as the Glauber formula, being a gener- 
alization of Glauber's result for h-A scattering, which has been derived, however, 
within nonrelativistic scattering theory. Eq. (78) is still very complicated, being 
a AB-fold integration, which cannot be reduced to a product of integrations as in 
the h-A case. Eqs. (73) and (78) are our basic formulas for the Monte Carlo treat- 
ment, which allow to generate appropriate collision sequences and thus to simulate 
nucleus-nucleus scattering. 
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